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Abstract 

We study the ground state <l/o = ^2 a x \X) of TV hard-core bosons on a finite lattice in configu- 
ration space, X — {xi, . . . ,xn}- All ax being positive, the ratios ax/^ay can be interpreted as 
probabilities P a (X). Let Eq denote the energy of the ground state and \dX\ the number of nearest- 
neighbour particl e-hol e pairs in the configuration X. We prove the concentration of P a onto X 
with \dX\ in a * / 1 -Eo |-neighbourhood of \Eo\, show that the average of ax over configurations with 
\dX\ — n increases exponentially with n, discuss fluctuations about this average, derive upper and 
lower bounds on Eo and give an argument for off-diagonal long-range order in the ground state. 
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1 Introduction 



Bosons on a lattice interacting via an infinite on-site repulsion (hard-core bosons) represent a system of 
double interest. This is the simplest example of an interacting Bose gas and, thus, the most promising 
candidate for a rigorous treatment of Bose-Einstein condensation of interacting particles. The model 
is also known to be equivalent to a system of h spins (Matsubara and Matsuda 1956) coupled via the 
X and Y and possibly the Z components of neighbouring spins and exposed to an external magnetic 
field in the Z direction. Ordering of the planar component of the spins is equivalent to Bose-Einstein 
condensation or the appearance of an off-diagonal long-range order (ODLRO, Yang 1962) in the system 
of bosons. Apart from some exceptions, such as bounds on the density of the condensate (Toth 1991) or 
the discussion of the model on complete graphs (Toth 1990, Penrose 1991), the most interesting rigorous 
results were formulated in spin terminology and obtained by the use of a particular symmetry, reflexion 
positivity (Dyson et al 1978, Kennedy et al 1988, Kubo and Kishi 1988). This introduced limitations on 
the value of the external field (zero field) and the lattice type (essentially hypercubic lattices). Expressed 
in terms of a boson gas, ODLRO was shown only at half filling on hypercubic lattices: in the ground 
state in dimensions > 2, and for positive temperatures above two dimensions. The proof of ordering does 
not offer much insight into the structure of the state. The only case in which details are known is the 
mathematically nice but physically trivial example of the complete graph (Toth 1990, Penrose 1991). 

In this paper we provide informations concerning the ground state \I/o> valid for very different lat- 
tices and arbitrary particle densities. The main result is the proof of a strong non-uniformity of in 
configuration space. It appears mathematically in the form of a large-deviation principle and makes the 
ground state resemble a thermal Gibbs distribution of a classical gas on the same lattice. It also leads 
to an approximate expression for and to an argument for Bose-Einstein condensation in the ground 
state in two and higher dimensions. 

Let L be an infinite lattice which, for the sake of simplicity, will be supposed to be regular with a 
constant coordination number (degree) k. Throughout the paper A denotes a finite connected part of L 
taken with periodic boundary conditions so as to keep the degree constant (not really essential). The 
Hamiltonian we are going to study in detail is 

Ho = -J2(Kby + b* y b x ) . (1) 

(xy) 

We write x,y, . . . for the vertices (sites) and (xy) for the edges (nearest-neighbour pairs); b* x and b x create, 
respectively, annihilate a hard-core boson at x. Boson operators at different sites commute with each 
other while 

b* x b x + b x b% = 1 (2) 

accounts for the hard-core condition. Correspondence with spin models is obtained by setting b x = S~ 
and b* — S x . The Hamiltonian conserves the number of bosons, 

and is also invariant under particle-hole transformation. We can, therefore, fix N so that p = A/|A| 
is between and \. (Here and below, if A is a finite set, \A\ stands for the number of elements.) Let 
X,Y, . . . denote A-point subsets of A, also called configurations. A convenient basis is formed by the 
states 

i*>=(n^)i°> ( 4 ) 

Vex / 

where |0) is the vacuum state. According to the Perron- Frobenius theorem, there is a single ground state 
whose all coefficients in the basis ([!]) are positive. Let a — (ax) denote the vector of these coefficients. 
Since the unique source of energy is hopping constrained by an on-site exclusion, we can expect that ax 
increases, at least 'on average', with an increasing number of nearest neighbour particle- hole pairs. More 
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precisely, let dX denote the set of nearest-neighbour particle-hole pairs for X (the 'boundary' of X) and 
Q, n the ensemble of configurations X with 'boundary length' \dX\ = n. It is easily seen that 

\dX\ = (X\H 2 \X) = (X\A\X) (5) 

where 

A = 2^[n x (l - n y ) + n v {l - n x )} (6) 
(xy) 

is the diagonal part of Hq, and |Q„| is the multiplicity of the eigenvalue n of A. We expect that the 
X-dependence of ax appears mainly through \dX\, and the average 



is some increasing function of n. We shall indeed find that a n rapidly grows with n - exponentially fast 
for n of the order of N and even faster for smaller n. To make a comparison with the free Bose gas, let 
us note that the ground state of the latter reads 



* free = |A|- \ L m , n«) m3: i o > w 

where a* is the ordinary boson creation operator at x. If we project out the part with no particle 
encounters, we find constant x \X), the sum going over the basis (^J). This uniform sum is qualitatively 
very different from what we actually find for the ground state, and the energy of this state will also be 
found to be extensively higher than that of the ground state. 
We start, in section 2, by studying the probability measure 

Pain) = . (9) 

Lm \ il m\a m 

We show that the mean value according to p a is \Eq\, the modulus of the ground state energy, while the 
mean square deviation D 2 is of the order of N, so that p a (n) is concentrated onto integers with \n — \Eq\\ 
of the order of \/N . In section 3 we evoke a large deviation principle for the Ising model on the same 
lattice, having the approximate form 

q{n) EE '""I ex e-^ 1 . (10) 



Here M = \dX\, the (arithmetic) mean of boundary lengths among TV-point configurations (and also 
the modulus of the energy of the projected free-boson ground state) and D 2 is the corresponding mean 
square deviation of \dX\. We derive a formula for D 2 as a function of p and k (M = kp(l — p)\A\ can 
be obtained trivially). In section 4 we use the results on p a and q to give an approximate expression 
for a n in terms of Eq, D a , M and D. We present an argument for the monotonicity of a n and deduce 
from it upper and lower bounds on D a /D, predicting a n — constant if \Eq\/M — 1 asymptotically 
(as |A| and N go to infinity). The variational treatment in section 5 shows that \Eq\/M > 1 in the 
thermodynamic limit. A notable exception, with |_Eo| = M = N(\A\ — N), is the complete graph where 
by permutational symmetry we actually have ax — constant. In section 5 we also present lower bounds 
on Eq for all densities. The physical consequences of our findings are resumed in section 6. Here we 
provide an argument for the existence of ODLRO in the two and higher dimensional models. It is based 
on a hypothesis which also permits to estimate the dimension dependence of the deviations of ax from 
the average a\gx\ and to conclude that these deviations are irrelevant above two dimensions. In section 7 
we briefly discuss extensions to Hamiltonians with more complicated interactions. 
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2 Non-uniformity of the ground state 



The results of this section are valid for any connected, not necessarily regular lattice with all coordination 
numbers < k. 

Let A = (Axy) denote the matrix of —Hq in the basis (|J). Then A is the adjacency matrix of a 
graph G whose vertices are the configurations, and two configurations X and Y form an edge of G if and 
only if their symmetric difference is an edge of A: X\YUY\X = {x, y} = (xy). We denote the largest 
eigenvalue of A by Ai and the corresponding eigenvector by a = (ax)- The ground state energy and wave 
function are related to them through Eq = —X% and = X) a x\X). The coefficients of the ground state 
wave vector a being positive, ax / ^2 ay can be interpreted as a probability P a (X) of X . As we shall see, 
P a (X) is very far from being a uniform distribution that one can find on complete graphs. 

A well-known fact concerning adjacency matrices is that (A n )xY provides the number of walks of 
length n in G between X and Y (Biggs 1974). Let 

W n {X) = J2( An )xY , (11) 

Y 

the number of walks of length n starting from X. First, we note that the expectation value of this number 
with respect to P a , 

(w n ) a = j2w n (x)a x /j2*x = x? = {w 1 ) n a • ( 12 ) 

x x 

This follows by taking the X component of the vector equation A n a = A™a, summing over X and dividing 
by ^2 ax- The above equality holds true for any simple graph: the probability measure P a 'sees' the 
graph as if it were regular with a degree Ai. 

Next, we make use of the fact that actually we are dealing with a sequence of graphs, G = Ga,jv, 
having the particular property that the typical degree is of order N while its change between neighbouring 
vertices is of order one: The degree of X is its number of neighbours, W\(X). There is a one-to-one 
correspondence between nearest neighbour particle-hole pairs if X is considered as a subset of A and 
neighbours of X as a vertex in G. Therefore, dX can be identified with the set of neighbours of X in 
G and we have W\{X) = \dX\. When passing from X to a Y S dX a neighbouring particle-hole pair 
is interchanged. For both the particle and the hole the number of neighbours of the opposite kind can 
change by at most k — 1, whence 

\Wi{X) - Wi{Y)\ < 2(k- 1) if YedX. (13) 

Now 



W 2 (X)= £ \8Y\ (14) 

YedX 

which together with (|l^) yields 

IW^X) 2 - W 2 (X)\ <2(k-l)\8X\ . (15) 
Taking the expectation value and using (|l^), 

0<D 2 a = (W?) a -(W 1 )l<2(k-l)X 1 . (16) 
By Chebyshev's inequality we then find that for any e > 

Pa (\Wi - Ai| > y/2(k - l)Ai/e) < E (17) 

or equivalently 

a x >(l-e)Y,ax- (18) 

|Wi(X)-A l |< 1 /2(fc-l)A l /e 



4 



Similarly to (14) 



w n {x)= E E • • E 



which, together with (fljf), yields 



n-l 



W n < J| (Wi + 2(k - 1)1) W n > Y[ (Wi - 2(fc - 1)0 



(19) 



(20) 



1=0 



1=0 



These bounds are nontrivial if n <C |A|, for example, for n fixed and |A| going to infinity. (In the opposite 
limit we have the stronger inequality W n (X) < (kN) n .) From (20) we get 



\W?-W n \ < J](2fc-2)' 

m— 1 

Taking the expectation value, 

n-l 

KWI%-A?| < E( 2fc - 2 ) r 



. l<h<---<l m <n-l 



(21) 



E 



if-i m \ (wr m ) a 



(22) 



(23) 



m=l \l<h<—<l m <n— 1 

Replacing n by n - m we obtain 

(w? _m )a < Xi~ m + o^-" 1 - 1 ) . 

Now Xi is of order |A|, so finally 

KW^a-X^^nin-lKk-m+OdM- 1 )}^ 1 ■ (24) 

Equations ([i6|)-(|l8|) and (|J) describe the concentration of the probability measure P a on configurations 
X whose boundary length \dX\ is in a \/iV-neighbourhood of \Eq\. We shall refer to this property as the 
non-uniformity of the ground state. Summing over X in Q„ , P a gives rise to the probability distribution 
p a , which is, hence, peaked about \Eq\. In section 4 we shall arrive at the same conclusion in a different 
way, by studying first the a priori distribution q(n) = |f2 n |/ ^2 m \Q m \- 

3 Large-deviation principle for the a priori distribution 

The precise computation of |f2„| is a difficult and unsolved combinatorial problem. The logarithm of 
this number is the entropy of the Ising model, for energy n, in a microcanonical ensemble with a fixed 
magnetization, ^2 xe \Cx = 2A^ — |A|. Indeed, if we put o~ x = 1 for x in X and a x = —1 elsewhere, 
we get an Ising configuration with a density p of + spins. The corresponding union of contours can be 
identified with OX whose total length \dX\ is the energy of the Ising configuration. For such systems 
there exists a strong version of the equivalence of ensembles which, applied to our case, states that the 
distribution q of |<9A| satisfies a large deviation principle whose rate function is, apart from a shift, minus 
the specific entropy s(e, p) of the Ising model (e = n/|A|), cf. Pfister 1991, Dobrushin and Shlosman 1994 
and Lewis et al 1994. The exact form of the entropy is unknown. To circumvent this problem, we use an 
approximate formula for the probability of having a boundary of length n, 



i r (n- M) 2 _ 
q(n) a Z- 1 exp{- 1 . ; } 



2D 2 



(25) 



Here Z oc y N is for normalization and 



M = \dX\ = kp(l - p)\A\ D 2 = {\dX\ - \8X\) 2 



(26) 
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cf. Eq. ([To]) in the Introduction. The Gaussian approximation is correct in a neighbourhood of the 
maximum of the specific entropy if |A|/D 2 is non- vanishing in the thermodynamic limit. This is what 
we are going to check by explicitly computing the mean square deviation of \dX\. We show that 

D 2 = [k - 2{2k - l)p(l - p)]kp(l - p)|A| = [k- 2(2k - l)p(l - p)]M . (27) 

We have found this expression first for d-dimensional hypercubic lattices, and have also checked it for 
the triangular, honeycomb and Kagome lattices. The apparent independence of the details of the lattice 
is somewhat surprising because our derivation below needs knowledge of the rather different local neigh- 
bourhoods up to next-nearest neighbours. The only common feature of all these lattices seems to be that 
all sites are symmetry- related and thus equivalent. This alone should therefore suffice to prove Eq. (p7j). 
Let us note that, with the microcanonical temperature defined as (9s/<9e) _1 , the maximum of the entropy 
corresponds to infinite temperature. Therefore, in a neighbourhood of the maximum we are at a safe 
distance to the ferro- and antiferromagnetic phase transitions, which could weaken our approximation 
through a breakdown (in finite volumes) of the concavity of the microcanonical entropy (Pleimling and 
Hiiller 2000). 

Equation (|2^) is obtained by filling the sites of A independently and with equal probability p. We 
expect smaller order corrections to appear if the computation is done with N fixed, N/\ A| = p (cf. Eq. ( p2[ ) 
below). So in this section X is a random subset of A whose probability to be selected is p\ x \ (1 — p)' A '~^T, 
and n x = n x (X) is a random variable taking the value 1 if x is in X and otherwise. Then all n x are 
independent and take 1 with probability p and with probability 1 — p. 

We define 

L = n x (! - n v) ( 28 ) 

y€dx 

where dx denotes the set of neighbours of x in A. The boundary length of X is obtained (cf. Eqs. (||) 
and (||)) as 

\dX\ = Y,UX). (29) 

xeA 

Thus the mean value of (|29|) is 

M = Y,U = kp(l-p)\A\ (30) 

x£A 

as claimed earlier. 

Let d(x, y) denote the graph distance of x and y in A, i.e. the length of the shortest walk between 
them. Since f x and f y are independent if d(x, y) > 2, we find 



D 2 = J2(f*-fM-fv)= E r ^y) ( 31 ) 

x,y x,y.d(x,y)<2 



r(x,y) = fxfy- fx ■ (32) 

The computation of the different terms is straightforward by observing that n x — n x . The contribution 
of the diagonal terms x = y is the same for any /c-regular lattice. Namely, 

Jl =kp(l-p) + k(k-l)p(l-p) 2 (33) 

r(x, x) = \A\r(x, x) = [k - (2k - l)p + kp 2 ]M . (34) 

The contribution of nearest neighbour pairs depends on the number of triangles containing a given edge. 
If there are I such triangles then 



fxfy = P 2 P) + [(k I? I } (1 - p) 2 ] (35) 
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r(x,y) =k\A\r(x,y)=p[£p- (2k - - p)]M . (36) 

x,y: d(x,y)=l 

If x and y are next-nearest neighbours to each other, they may have m common nearest neighbours. 
Then 



fxfy = P 2 [m(l-p) + (k 2 - m )(l-p) 2 ] 
r(x,y) = mp 3 (l-p). (37) 

In d-dimensional hypercubic lattices (k = 2d) there are next-nearest neighbour pairs with m = 1 and 
m = 2. Their contribution to D 2 is 



r{x,y) = k\A\r{x,y) m=1 +4( d )\Mr(x,y) 



x,y: d(x,y)=2 



m=2 



2j 

p 2 M + (k-2)p 2 M = (k~l)p 2 M . (38) 



For the triangular lattice (k = 6) 



r(x,y)=k\A\[r(x,y) m=1 +r(x,y) m=2 ] = (k-3)p 2 M . (39) 

x,y: d(x,y)=2 

In the honeycomb lattice (k = 3) each site has six next- nearest neighbours, all of the type m = 1. So 

]T r(x,y) = 2k\A\r(x,y) m=1 = (k-l)p 2 M . (40) 

x,y: d(x,y)=2 

In the Kagome lattice (fc = 4) there are eight next-nearest neighbours with m = 1: 

^ r(x,y) = 2k\A\r(x,y) m=1 = (k-2)p 2 M . (41) 

x,y. d(x,y)=2 

Finally, we obtain D 2 by adding (|3j) and (|36| ) with i = and (|3^) for hypercubic lattices, (|34| ) and 
( ^6| ) with £ = 2 and ( ^9|) for the triangular lattice, and ((3^) with £ = and (^) for the honeycomb 
lattice and @, © with ^ = 1 and © for the Kagome lattice. All yield @. 

If the random variables were independent, the mean square deviation of their sum would be given by 
(|34|). For any p <\ this is larger than the actual value (|27|), so on average the f x are negatively correlated. 
Because I < k — 1, the nearest neighbour correlation r(x,y) < for all p < |, cf. Eq. (|3^). On the other 
hand, according to (|3^), for next-nearest neighbours r(x,y) is always positive. Apparently, the induced 
probability distribution for the f x could be approximated by a Boltzmann-Gibbs distribution written 
with an antiferromagnetic nearest-neighbour Hamiltonian for the centralized variables f x — fx and with 
a temperature and an external field chosen so as to fit the computed average (^0|) and nearest-neighbour 
correlations ( |36| ) . If the effective j3 and field are below their respective critical values we get exponentially 
decaying antiferromagnetic correlations, and the approximation is qualitatively correct. Adding next- 
nearest- neighbour interactions one could fit the computed next-nearest-neighbour correlations, and so 



4 The form of the averaged wave vector 

Using the large-deviation result for q(n) we can obtain more precise information on p a (n) and a n . Apart 
from normalization p a (n) is obtained by multiplying q(n) with a n . Since, as we shall see in the next 
section, Ai > M and the difference is of the order of N, a n has to increase exponentially fast - at least 
in a neighbourhood of M - so as to shift the expectation value M of q to the expectation value Ai of 
p a . As a result, we obtain the approximate expression p a (n) ~ exp— (n — Xi) 2 /2D 2 , consistent with our 
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findings in section 2. This is a second, independent, argument for the non-uniformity of the ground state, 
assuming the form of another large-deviation principle for \dX\. To get a n we equate this form of p a (n) 
with the other one, (^|), where |fi„| is estimated through (|l0|). This gives 

a ~ exp ( {n ~ M)2 - t^ll) K exp < g « = D ^ ± ^ = ^ (42) 



For the moment, we have only trivial estimates for Ai (> M) and D a (cf. Eq. (|lq)). More precise 
bounds on D a can be obtained from the monotonic increase of a n . 

The argument telling that a n should increase for all allowed n is as follows. From the spectral 
decomposition of A m we obtain (supposing J2 a x = 1) 

a x a Y = lim X^ m (A m ) XY (43) 

m — >-oo 

if A is non-bipartite, and 

a x a Y = lim I[Ar" l (A m ) xy + X^' 1 [A m+1 ) XY ] = lim ^(^) xy (44) 

if A is bipartite; r m = 2m or 2m+l if the graph distance da(X, Y) of A and Y is even or odd, respectively. 
Summing over Y and averaging over X in f2 n we obtain, for example, for a non-bipartite A, 

a n = (£ a,)" 1 Jm A r ™-^ ^ W m (X) . (45) 

Inspecting Eq. (^9|) one can conclude that for any m the average of W m (X) over Q„ increases with n 
and, hence, a„ > a n ' if n > n'. 

Let n m ; n and n max denote the smallest and largest allowed values of \dX\, respectively. Since n m j n = 
o(N), monotonicity of a n implies through ( [42| ) 

nmax ~ Al < 1% < (46) 

up to an error of o(l). Therefore, the coefficient of n in the exponent of a n is non-negative. Clearly 
Xi/M = 1 would imply D a /D = 1 and a n = constant. In the next section we show that for finite 
dimensional lattices X\/M > 1 asymptotically and, hence, the two bounds in ( f46| ) form an interval 
around 1. On bipartite lattices n max = kN, so at half filling n max = 2M and 1 is in the centre of the 
interval. It is possible that D a = D and the exponent of a n is linear in n. However, even if D a =/= D, 
the quadratic term in the exponent is of the same order as or even smaller (if n = o(N)) than the linear 
term. 

More information can be extracted from Eqs. (|43| ) and (fl4|) if we combine them with (|l8|). If m is 
large enough, 

l - 2 ^W^X) E (A- )XY <1, (47) 



||ar|-A 1 |<^2(fe-i)A 1 / £ 



that is, for any X and any sufficiently large m an overwhelming majority of walks of length m starting 
from X on the graph G end up in a vertex whose degree is Ai + 0(>/~N). In the expression 

W m (X) = (a x ar)X? + o(X?) (m oo) , (48) 

obtained from p^), A" 1 accounts for this 'long-time' behaviour. Let us note that ( f47j ) cannot be un- 
derstood by imagining the walks on G as realizations of some simple random process. Because G is 
nonregular (unless A is a complete graph) , no Markov process can assign equal probabilities to all walks 
of equal length. As an example, for a locally unbiased random motion (which chooses among neighbours 
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with equal probability) the most probable individual walks of a given length are the 'descending' ones, 
those going towards vertices of a lower degree, and the less probable the 'ascending' ones. Because of 
the form ( |25| ) of the distribution of degrees, an X with \dX\ < M has typically more ascending than 
descending neighbours and vice versa for \dX\ > M. Therefore, the most probable degrees of end-points 
of very long walks would be nevertheless close to M (but not to Ai). 

Although the monotonicity of a„ holds for all n between n m ; n and n max , the validity of ( (42| ) is limited 
to a neighbourhood of the interval \M , Ai]. In particular, the growth of a n for n = o(N) is much steeper 
than the exponential predicted by (|42|). Indeed, let us write the eigenvalue equation in the form 

i^i £ «* = \!k\ ax ■ (49) 

1 1 YedX 1 1 

For a d-dimensional lattice n m i n = 0(N^~). If we choose \dX\ — n to be of this order, we find that 
the average of ay over dX yields N x / d times ax! Since q(n) is also rapidly increasing here, the average 
over dX is dominated by Y with n + 2 < \dY\ < n + 2k — 2. From this and the monotonicity of a n we 
conclude that a n+ 2k-2 / ct n is at least of order N 1 ^. 

5 Bounds on the ground state energy 
5.1 Variational upper bounds 

Variational estimates of the ground state energy are of the form 

E < (i>\H\^)/(i>\i>) . (50) 

A trivial choice is 

</> = £|AQ (51) 
with the summation going over all iV-point subsets of A. It yields Eq < —M where 

2\EG\ 2\EK\^ZI 



M=\dX\ = -^= -^± = k (l-p)N + 0(l) , (52) 

cf. Eq. (|30|). Here \EG\ and \VG\ denote the number of edges and vertices of G, respectively, and \EA\ 
the number of edges of A. It is not a priori obvious that this bound can be improved in the order of the 
volume, and it is important to know that it really can. In the opposite case, if —kp{l — p) were the true 
ground state energy per site then, as in the complete graph, the Hamiltonian (Q) would have a product 
ground state in an infinite volume, 

^n^K^i+^K^)) , (53) 

X 

i.e. no local perturbation could decrease the energy of 'J. (Because of the product structure, if the energy 
could be decreased locally, the specific energy —kp(l — p) of ^ could also be decreased.) Since ^ shows 
ODLRO with the value of the order parameter at its theoretical maximum (p(l — p), cf. Appendix), by 
proving it cannot be a ground state we exclude a trivial scenario for Bose-Einstein condensation. 
To prove that |-Eo| — M is of order N we apply trial functions of the form 

i> v =Y, v x\X) v x =v(\dX\) (54) 

i.e. vx depending on \dX\ only. The variational bound ( p0| ) reads Ai > B(v) where 

B[v) = = EJWK\ ■ (55) 



9 



Here Ni(X) is the number of those neighbours of X having a boundary length \dX\ + 2i. In Eq. fl5q ) we 
have used ( |l3| ) and the fact that the parity of \dX\ is the same for all X: even if N is even and that of 
k if N is odd. 

The form ( [f2| ) of a n suggests that the best choice for v(n) would be e xn+yn with x of order 1 and 
y of order However, we do not expect the quadratic term to yield a significant improvement, and 
choose a simple exponential v(n) = e xn with < x < x max where 

"•max - M , . 

(56) 



2D 2 

Then by making use of the large-deviation principle for \dX\ we find 

B(v) = (n) x (e 2xi ) [(nh] (57) 

where 

w. - |fSw </<■».. - £ /o^j E »>m ■ (58) 

^ 1 i=-k+i 1 1 xea„ 

Now, since (n)o = M and d(n) x /dx \ x =o = 2-D 2 , we find 

(ra)* = M (1 + 2xD 2 /M + 0{x 2 )) . (59) 

We note that the approximate form ( |To| ) would yield the same result. 

Next, we turn to the second term of B{v). Obviously (e 2xl ) n > e - 2 ( fc - 1 ) ;c . This yields B(v)/M > 1 
for p near but not near |. However, we can use Jensen's inequality {e 2xl ) n > e 2x ^ n together with the 
fact that (i)[Ml goes to zero as the volume increases. We then conclude that asymptotically (i) [{ n } x ] — 
— constant x x + higher order terms, so that the average of the exponential in Eq. ( |57|) equals 1 — 0(x 2 ). 
Thus, for x small enough we indeed find B(v)/M > 1. 

With somewhat more effort one can actually compute a lower bound on B(v). We do not make it 
here, only notice that in optimizing such a bound according to x (and also in optimizing a lower bound 
on Eq, see below) the knowledge of n max is necessary. On bipartite lattices n maK — kp\A\ for all p < g, 
and it is an easy graphical exercise to see that the same equality holds for p < | on the triangular and 
Kagome lattices. However, for both lattices n max is constant between the densities ^ and |: 2|A| for the 
triangular and ||A| for the Kagome lattice. This can be seen from the following argument. In general, 
— n max is the ground state energy of the antiferromagnetic Ising model under the restriction that the 
magnetization is fixed, X^eA 17 ^ = i^P ~ 1)1^1- However, we do not need to deal with the restriction. In 
both cases the (unrestricted) ground state is known to be highly degenerate. Among the exponentially 
large number of ground state configurations there are non-magnetized ones, corresponding to p = g, 
others with a concentration of up-spins p = i, and between these two limits p can vary by steps of 1 / 1 A | . 
The rule is to flip zero-energy spins one by one. The common energy of all these configurations is easy to 
compute from the fact that in each triangle there is precisely one unsatisfied bond. This fixes the value 
of n max as given above. 

We also note that the optimal x (which maximizes B(v)) is small: of order 0.1 or smaller. We shall 
evoke this fact in the discussion of section 6. 

5.2 Lower bounds 

In their paper Dyson, Lieb and Simon (DLS) gave a lower bound on the ground state energy of the spin-^ 
XY-model (Dyson et al 1978, Theorem C.l). Because it corresponds to an upper bound on the norm 
of the Hamiltonian (|l|) in Fock space, it is automatically valid for the hard-core boson gas (|]) at any 
density. It reads 

|JSb| < W VHkTY) if* is even 

1 1 ~ 4 1 1 \ k + 1 if k is odd . v ' 
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The above bound is to be compared with the trivial bound \Eq\ < n max obtained by writing Eq. ( fl9| ) for 
an X = Xq which maximizes ax- We then see that (B0j) is nontrivial for the d-dimensional hypercubic 



lattice if p > | J 1 + ^ , for the triangular lattice if p > , for the honeycomb lattice if p > | and for 

the Kagome lattice if p > \\J~^- 

The bound ( |60| ) is the best at p = \ and can be improved for lower densities. 
If eo is the ground state energy per site then 

M - (*o|A x |* ) = (*o|i E ^ b v + WI*o) (61) 

yGdx 

for any x in A. A x preserves the number of bosons in = {x} U dx, therefore it commutes with P x ,j, 
the orthogonal projection to the subspace Tt x j of states with j particles in A x , where j = 0, 1, . . . , k + 1. 
We can write 

(»,ia.i»o) = E ^wr ^ 01 ^ 1 ^ ■ (62) 

Since (^ol-fx.j ' ^rjl^oV^ol-Pzj'l^o) is a normalized positive linear functional, 

(*o|A x |* ) < E Amax (-?X*ol P -jl*o) ' ( 63 ) 
j 

where A max (j) is the maximum eigenvalue of A x restricted to H x ,j- Now A max (j) was computed by DLS: 

ArnaxO) = \Vj(k+l~]) ■ (64) 

The DLS-bound ( |60| ) corresponds to the maximum of (|64|). Let /i(t) be obtained by linear interpolation 
through the points (j, A max (j)): 

h(t) ^l{([t\+l-t)V[t\(k + l-[t\) + (t-[t\W([t\+l)(k-[t\)} . (65) 

Since h(t) is a concave function, and the right member of ( |63| ) can be written as J h(t)dis(t) with a 
probability measure i/ concentrated on the integers from to k + 1, by Jensen's inequality 

\e \<h(^J tdv(t)j = h ^E^ol-Pa^o) I • (66) 

The expectation value in the argument of h can be evaluated and yields (k + l)p. Indeed, using the 
translation invariance of 'J'oj after a simple algebra we find 

Ej<*o|p«j|*o) = E a *mEl* nA *i ■ ( 67 ) 

Let \x and xa x denote the characteristic functions of X and A x , respectively. Then 

\xnA x \ = ^2xx(v)xaM (68) 

0jiT.d tlioi'oforo 

Ei xnA -i = E^(y)E^^) = i x i-i A -i = 7V ( fc + 1 ) ■ ( 69 ) 
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Figure 1 : Approximations of the ground state energy per site as a function of the density in the square 
lattice. 

Hence, we obtain 

M < mm{h((k + l)p), kp} < i(fc + l)Vp(l-p) (70) 

where the last member results from majorizing h(t) by \\/t(k + 1 — t). Taking the minimum in Eq. ( f70| ) 
is not superfluous. For p < we have h((k + l)p) = jvk(k + l)p. Since n max /|A| = kp is smaller 
for all k, the trivial bound is better for small densities. Actually the trivial bound provides the right 
asymptotics at p = because eo = — kp + 0{p 2 ) near p — 0, as one can see by comparing the trivial 
upper and lower bounds, —M/\A\ and — n max /|A|. 

Numerical works on the spin-^ XY model in the square lattice yield an approximate formula, 

eo (p) = -1.09766 + 4.835(0.5 - pf - 1.99232(0.5 - p) 4 + 0.85952(0.5 - pf . (71) 

The first two terms come respectively from a Monte Carlo (Zhang and Runge 1992) and a finite-size 
scaling (Hamer et al 1999) study, to which we have added the fourth and the sixth order terms to obtain 
ep(0) = and ed(0) = —4. In Figure 1 we have plotted the trivial upper bound — 4/o(l — p), the formula 
(fry) and the lower bound — min{h((k + l)p), kp}. We have also shown numerical points for a 5 x 5 lattice 
from Table 1 of Hamer et al 1999. They nicely follow the curve (|7l|). 

Using the formula (|7l|) we can compute the lower and upper bounds ( |46| ) on D^/D 2 . Both of them 
increase monotonically as p varies between and 0.5: the lower bound from ~ 0.66 to 0.90234 and the 
upper bound from 1 to 1.09766, their difference decreasing monotonically at the same time. 



12 



6 Off-diagonal order and the final form of the ground state 



Off-diagonal long-range order in the ground state is characterized by (cf. Appendix) 



CJ 



= lim ^j2E a * E E 



axu{y}\{x} 



> 



(72) 



A" 



xex y £x 



where lim means the thermodynamic limit (|A| and N going to infinity and A/|A| going to p). In ([72| 
^a 2 x — 1 is supposed. 

Let us see, what would we obtain if ax did not fluctuate within f2„, i.e. if ax = Q-\dX\- This and (|42| 
would imply another large-deviation principle in the approximate form 



exp 



(n-X 2 ) 2 



2D 1 



with 



A, 



Ai + D 



2 Ai -M 
' a 2D 2 -D 2 

The mean value A2 is larger than Ai because 

Dl 



n 2 - 



D Dl 



2D 2 - D 2 n 



Ai 

„ < — < mm ■ 
D 2 ~ M 



1 



fc+1 



1 



P 2k^p{l-p) 



< 2 



see Eq. M$). As 



\\d(XU{y}\{x})\-\dX\\<2k 
with a XU { y }\{x} > a\ dX \-2k we would find 



w o > p(l — p) lime 



2fc(A 1 —hi) 
— r>^ 



> 



(73) 
(74) 

(75) 
(76) 

(77) 



The argument above, together with the conclusion (|77|) may be right in high enough dimensions but it is 
known to be incorrect in one dimension. Fluctuations have to be taken into account. 
A way to prove (f72|) would be to show 



a xu{ y }\{x} > ca-x 



(78) 



for all A, x and y with a c > independent of A and N . If (|78|) did hold true, we would find luq > 
cp(l — p) > 0. In the complete graph ([78]) is verified with c = 1 and equality sign. However, we know 
already from the considerations following Eq. (^9|) that in finite dimensional lattices (|7^) cannot hold 
indeed for all X: e.g., in the case when \dX\ cx A 1_ 3 ; ( |78|) can be satisfied only with c cx N~i . Also, 
we should not be able to prove ([78]) in one dimension. Below we present an argument which takes into 
account fluctuations and distinguishes between one, two and higher dimensions. 

First, we note that we need (78) only for A in a subset S(A, N) such that the sum of a 2 x over S is non- 
vanishing in the thermodynamic limit. Let us introduce a function R{X) by setting ax = &\dx\ eR< " X ^ ■ 
Then 

£ e«W = |O n | < £ e 2R ^ < \n n \ 2 , (79) 

xen n xen n 



so we can define an e n such that < e„ < 1 and 

e 2Ri - x) 

xen n 



I fin 



(80) 
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Using this and Eq. we find 



E 



' A" 



exp 



(n - M n f 
2DZ 



with 



M n = Ai + — i 1 /" )£> ° „ (Ai - M) 
2I? 2 - (1 - cO-DT 



D 2 = 



D 2 D„ 



" 2£ 2 - (1 - e„),D 2 



(81) 



(82) 



Because of the possible n-dependence of M n and D n , Eq. ( plj ) may not describe a large deviation principle. 
However, M„ and D n satisfy the inequalities Ai < M n < A2 and D a /^/2 < D n < D a i and, hence, for 
any e > 

~ • (83) 



lim 



E 



' A 



(l-e)\ 1 <\dX\<(l+e)\ 2 



For S we can take the set of summation with any e > 0. 

Next, we observe that by using Eqs. (B3I) and (E3) we can write 



ax ,. 
= hm 

ay 



m^TSo W m (F) 



lim 



A 1 W m (X) + VF m+1 (X) 



(84) 



for non-bipartite or bipartite lattices, respectively. The convergence in (|84j) is usually slow. However, we 
expect that if \dX\ and |<9y| do not differ too much then, at least up to the order of magnitude, ax jay 
can be approximated by the ratio of the number of walks whose length is the distance between X and 
Y. Thus, we conjecture that for \dX\ = Cl N and \ \dY\ - \dX\ \ < c 2 



, a x 
In 



In 



W dGix ,Y)(X) 



ay W dGi x,Y)(Y) 



< c 3 



(85) 



where C3 may depend on the constants c\ and C2 but not on the size of the system. In what follows, we 
examine the consequences of this hypothesis. 

If the aspect ratios of A are kept bounded, its diameter L is of order |A| ^ and the diameter of G is 
of order L d+1 . Over such a distance ( |S5|) allows a gigantic change of ax- However, the distance to be 
considered for (|78| ) is much smaller. If the symmetric difference of X and Y is {x, y}, it is easily seen 
that dc{X, Y) = d(x, y) < L. Let \dX\ = n cx N. Above one dimension we apply the inequalities ( f20|) to 
obtain 



W dG(X ,Y)(X) 



W, 



< 



1 + 2(k - l)L/n 
1- [2(fc- l)L + 2]/n 



(86) 



'd G (X.Y)\ 

In the thermodynamic limit the upper bound remains finite in two dimensions and converges to 1 in 
higher dimensions, so we get the necessary estimate (f78|). In one dimension there exists no convergent 
upper bound. Thus, we find ODLRO in the ground state in two and higher dimensions but not in one 
dimension. 

The difference between 2 and 2 + e dimensions in the convergence of the right-hand side of Eq. ( |S6"|) 
indicates that d = 2 is the critical dimension above which the deviations of ax from the average a|gx| 
are irrelevant. This remark can be made quantitative if we observe an interesting consequence of the 
hypothesis ( |S5| ) on R(X). If \dX\ cx N then inside a sphere of radius cx VN, centred at X, for all Y such 
that ||<9y| — \dX\ \ is bounded by some constant of the order of unity \R(Y) — R(X)\ also remains below 
another constant of order 1. Comparing the scale ~ y/~N with the diameter ~ N 1+1 / d of G we conclude 
that 

\R{X)\ < mnxR(Y) - minify) < const x N^ + ^ . (87) 

The maximum and the minimum are taken over f2|gx|- The first inequality holds because for any n the 
average of R(X) over fi„ is negative, but R(X) cannot be negative for all X. 
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According to Equation (|8 

a x = exp[F(\dX\) + R(X)] = exp F(\dX\) + 

Therefore, in one dimension R(X) can dominate F(\dX\), in two dimensions it can vary on the same 
scale ~ N, and above two dimensions the variation of R(X) is negligible compared with that of F(\dX\). 
Apart from an additive constant necessary for normalization, the approximate form of F(n) = lna„ can 
be read off from Eq. ( ji^ ) . 

The upper bound in the right-hand member of Eq. ( |87| ) could, of course, overestimate the order of 
magnitude of the maximum of |i?(A)|. However, because it seems to correctly distinguish between one, 
two and higher dimensions, it is probably sharp. This implies also that on the scale ~ \/N the function 
R(X) can be considered as describing a ballistic motion over G rather than a random fluctuation about 
a negative average, which would lead to a smaller bound. Accepting ballisticity on smaller scales as well, 
one can find the order of magnitude of local variations of R(X). We distinguish between two cases. 

(a) If A is bipartite with sublattices Vi, V2 then G is also bipartite with sublattices Ti = {X : 
\X n Vi \ is even} and r 2 = {X : \X n V\ | is odd}. Along each path in G the subsequent terms belong to 
alternating sublattices. If d 2 X denotes the set of second neighbours of X then 

\R{X) - R(Y)\ < const x N~ 1/2 if Y e d 2 X and \dY\ = \dX\ . (89) 

(b) If A is not bipartite, neither is G, and no systematic compensation of terms of different signs is 
possible. In this case 

\R(X) - R(Y)\ < const x N~ 1/2 if Y e dX and \dY\ = \dX\ . (90) 



The result on ODLRO could be obtained from (|89|), (J90j) as well. 

As the two-dimensional case is marginal, some subtle logarithmic corrections, that we do not see, may 
modify our conclusions. 

Above two dimensions we also find, as a further consequence of Eq. (p7|) , that e n identically vanishes 
in the thermodynamic limit. Thus, M n /X2 — > 1 and D n /D a 2 — > 1, Eq. (pll) asymptotically coincides with 



Eq. (|T3h, Eq. (|77|) provides a valid lower bound on the order parameter, and Eq. (83) can be replaced by 



lim Yl a x = 1 ( 91 ) 

(l-e)A 2 <|0X|<(l+£)A 2 

for any e > 0. This property makes a 2 x similar to a thermal Gibbs state for a classical lattice gas. 
The classical Hamiltonian (with 1/fcgT incorporated) is H e g(X) = —2F(\dX\) — 2R(X), and there is 
asymptotic equivalence between the canonical distribution and a microcanonical one concentrated on 
X's with H c ff(X) in an o(N) neighbourhood of the canonical expectation value of H c g. This holds in 
spite of H c g containing long-range interactions. Indeed, — F is a sum of nearest neighbour repulsive 
two-body interactions, / x , and, if D a ^ D, properly normalized long-range four-body interactions, f x fy, 
see Eqs. (p8|), ( p9| ) and (^2|), the whole yielding a global repulsion which tends to maximize \dX\. Now 
R may also include up to N body interactions but because R(X) = o(N) and because a 2 x corresponds to 
a high-temperature Gibbs state, it can be neglected: The fact that R may partly split the ground state 
degeneracy of —F is irrelevant at high temperatures. 

We have no rigorous proof that the limiting effective Gibbs state is a pure (in particular, a high-tem- 
perature) state. For this we should show that for all densities and all relevant n, 2F(n)/n < K c , the 
critical coupling of the Ising model on the same lattice. Even if we knew Ai exactly, our bounds ( fl6| ) 
on D a /D are not good enough to obtain this information. For example, in the case of the square lattice 
we can use the numerical fit (|7l]). If we assume that D a > D, we get 2F(n)/n < 2(Ai — M)/D 2 where 
the upper bound increases with p and its maximum at p = 0.5 is 0.3906 < K c = 0.4407. On the other 
hand, if we use the lower bound for D a , at half filling we find for n > \\ that 2F(n)/n < 0.628 > K c . 
An indication to the high-temperature character can be found in the variational estimates, cf. the final 



15 



remark of section 5.1. In the light of Eq. (p8|), trial functions of the form ( p4| ) acquire a particular 
importance: For d > 2 one could find with them the exact ground state energy per site, and the optimal 
exponential ansatz may not be far from the true ground state. The variable x appearing in the formula 
( p7| ) corresponds to ^/3J of the Ising model. The estimates we have done predict a value well below 
the critical one in any dimension: [3J < 0.2 for all densities at k = 3, and decreases with increasing fc, 
reaching in the complete graph. 

Now the fact that a\ is equivalent to a high-temperature Gibbs state means that there is no classical 
(diagonal) order in the ground state, coexisting with the purely quantum-mechanical (off-diagonal) one. 
Diagonal order is characterized by an order parameter operator O which is diagonal in the basis (Q) 
and for which the ground state expectation value | A| — 2 (\E r |C , *0|\E f o) = | W | 2 J2 a x\0{X)\ 2 has a non- 
vanishing thermodynamic limit. A typical example is O q = ^2 X e % ^^n Xl associated with a periodic order 
with a wave vector q. Due to the correspondence with a Gibbs state, diagonal order in the ground state 
of the Bose gas is equivalent to low-temperature order in the associated classical lattice gas. The ground 
state of — F is exponentially degenerate for p < 1/2. However, it is only two-fold degenerate on bipartite 
lattices at p — 1/2, and the degeneracy is not split by the translation invariant R(X). Apparently, there 
is no qualitative argument against a crystalline order at half filling in the ground state of the d > 2 
dimensional purely hard-core Bose gas. The nearest neighbour repulsion in — F is simply not strong 
enough, reflecting the fact that a purely on-site repulsion in H , even though infinite, cannot induce such 
an order. 



7 Extensions and concluding remarks 

The nonuniformity of the ground state we have found in a purely hard-core gas seems to be a general 
property of finite-range Hamiltonians. The extension to more complicated hard-core Hamiltonians is 
immediate. Let 

H = -J2( b x b y + K b ^ + J2 v * n x- J J2 ( nx ~l) ( ny ~l) +K ^ n ^ 

(xy) (xy) V J V J 

= H a + V+^(2A-\EA\) + K (92) 

where A is the diagonal part of Hq, cf. Eq. (||), and K may contain further interactions of a bounded 
range. Because the terms additional to Hq are diagonal in the basis (||), the Perron- Frobenius theorem 
still applies and the ground state wave vector ax > 0. If A is the matrix of — H then keeping the 
definition ([ll]) we still have (|T2| ), and the inequalities of section 2 also hold with suitable modifications. 
In particular, in (jl^) and ( fl5| ) 2(fc — 1) is to be replaced by k — |2 — J\(k — 1) + SV + SK where 
SV = max( l;/ ) \V X ~ V y \ and SK = m&nx.Yedx \K(Y) — K(X)\ are of the order of unity, and in ( [l6| ) 
we have n(\dX\) a instead of 2(fc — l)Ai- Now 

W X {X) - A xx + \dX\ = (l - \dX\ - V(X) - K(X) + ^\EA\ (93) 

may not be positive and (Wi) a = Ai = — Eq may not be of order N. Therefore, there may not be a large 
deviation principle for the pair (W%, a). However, if we replace H by H c = H + cl where / is the identity 
operator, the ground state vector a will not change. The corresponding W\ is Wf (X) — W\(X) — c, and 
D^{Wi) is also independent of c. Let us choose c = (Axx)a, then 

A!+c=A^ = {Wf) a = {\dX\) a (94) 

and, thus, we have a large deviation principle for [W^a) in the same form as for the purely hard-core 
interaction. 

What really counts for the non-uniformity of a is not the large deviation principle for ( W\ , a) but 
the variation of Wi(A) over the set of A-point configurations. The ground state is nonuniform if 
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[maxWi(X) — min W\ (X)]/VN goes to infinity with N. As an example, let us consider the nearest 
neighbour anisotropic Heisenberg model (V = K = 0). At J = 2, when it is isotropic, for any fixed N we 
have Wi(X) = |i?A|/2 and hence Ai = \EA\/2 and ax is constant. Away from the isotropy point AWi 
is of order N and a is nonuniform. 

If il n — {X : W\{X) = w n }, where wi < W2 < ■ ■ ■ are the possible values of W\, one can write down 
an approximate formula similar to Eq. (|l0|), namely q(n) ~ exp[— (w n — W\) 2 /2D 2 (W\)](w n+ i — w n ). 
Interpreting a n as the average of ax over the newly defined Q n , we have also the analogue of Eq. ([42|). 

As earlier, we can write ax = exp[F(Wi(X)) + R(X)], and Eq. ( frit ) remains valid. The approximate 
form of F is 

F{Wl{X)) M 2/J(W/ 1 ) 2 2J a (^) 2 + ' (95) 

One can replace the infinite on-site repulsion by a finite one provided that stability is assured, H > 
— const x N. The dimension of the A-particle subspace will increase to (' A ' 4 ^ 1 )- A configuration X 
becomes a list of N elements containing vertices of A with possible repetitions that we can write as 
X = {cc nx } xe A- The neighbourhood relation remains the same: Y € dX if Y can be obtained from X 
by moving a single particle along an edge of A, say, from x into y. The corresponding matrix element 
of — H is then Axy — \J n x (X)n y (Y) . Now Y^y^x Ax"Y can vary between ky/N and \/2kN when \dX\ 
varies between k and kN. However, Axx can be of order N 2 and its change between neighbours of 
order N . Therefore, we miss a pointwise bound like ( |i"3| ) and cannot repeat the argument of section 2. 
Nevertheless, —Wi(X) = ^y(y|i?|X) is still the energy of a classical lattice gas in which the equivalence 
of the canonical and microcanonical ensembles implies the analogue of Eq. ( JTc| ) and, via the argument of 
section 4, the nonuniformity of the ground state. 

The real novelty introduced by the additional interactions is the possible coexistence of classical long- 
range order with the off-diagonal one. A thorough discussion of this question, concluding negatively, can 
be found already in Penrose and Onsager's 1956 paper. Nevertheless, a controversy has remained until 
recently, when reliable quantum Monte Carlo computations on square lattice models have shown such a 
coexistence (Batrouni et al 1995, Scalettar et al 1995). We hope that the ideas developed in this paper 
can contribute in the future to further elucidate this interesting problem. 

This work was supported by the Hungarian Scientific Research Fund (OTKA) under Grant No. T 
30543. 



Appendix. A digression on ODLRO 

Following the standard definition (Penrose and Onsager 1956, PO), we give below the expression of the 
ODLRO parameter in terms of {ax), valid for any Bose gas whose Hamiltonian contains a hard-core 
repulsion and conserves the number of particles. The density operator in the ground state ^0 (supposed 
to be normalized) is the orthogonal projection |\E'o)( v I / o|- The density matrix in the basis ([|) is therefore 
(ax ay)- From here the one-particle reduced density matrix a = (<J xy ) is obtained by taking a partial 
trace over N — 1 particle positions, 

°xy =/ j ax'u{x}0-X'u{y} ■ (96) 

X' 

In the above sum X' runs over the N — 1-point subsets of A not containing x and y. The matrix a is real 
symmetric, positive semidefinite (in fact, positive definite for N > 2), its trace is N and all its elements 
are positive. A possible choice for the order parameter, according to PO, is 

w = |Ap 2 °*v ■ ( 97 ) 

x,y£A 

Another definition of PO for the order parameter is 1/|A| times the largest eigenvalue of a. If A is a shift- 
invariant set (which supposes periodic boundary conditions) , the ground state is translation invariant and 
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&xy — @ x — 

y . Applying the Perron- Frobenius theorem to a we find that its largest eigenvalue is ^2 &x-y — 
|A|w, so that the above two definitions coincide. The order parameter u) also has an interpretation as the 
density of the Bose-Einstein condensate. This holds because the operator associated with the density of 
the condensate is 

lAI-^U^o - |A|- 2 b >y ( 98 ) 



and because 

Inserting (Q) into ( p7| ) we obtain 

|A|- ~ |A| 



* |&;M*o> =<t*v • (99) 



C =LJ ' + m2 J = JX\2 a * axu{y}\{x} ■ (100) 

X xeXygX 



The maximum of u>' under the condition of normalization is p(l — p), attained with ax constant. It is 
realized by Hq on complete graphs and by the ferromagnetic Heisenberg model on arbitrary graphs. Our 
ui' is one fourth of the order parameter usually used in the spin-^ XY model, cf. Fujiki and Betts 1986, 
Hamer et al 1999. There is ODLRO for a density p if u> or, equivalently, lj' has a nonvanishing limit 
uiq > as N and |A| go to infinity while iV/|A| goes to p. When this occurs, conservation of the particle 
number is spontaneously broken and there appears (at least) a one-parameter continuum of ground states 
in infinite volume. The order parameter operator, is nondiagonal in the natural basis (§) (whence 

the name ODLRO). The ground states 7 , < a < 2ir, are symmetry-related. For any local operator B 



.(B) = jo (l[e- ian *Bl[e ian A 

\ X X / 



(101) 



In particular, 

la (b* x ) = e~ ta MK) * . (102) 
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